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ABSTRACT 

The spin of a neutron star at birth may be impacted by the asymmetric character 
of the explosion of its massive progenitor. During the first second after bounce, the 
spiral mode of the Standing Accretion Shock Instability (SASI) is able to redistribute 
angular momentum and spin-up a neutron star born from a non-rotating progenitor. 
Our aim is to assess the robustness of this process. We perform 2D numerical simu¬ 
lations of a simplified setup in cylindrical geometry to investigate the timescale over 
which the dynamics is dominated by a spiral or a sloshing mode. We observe that 
the spiral mode prevails only if the ratio of the initial shock radius to the neutron 
star radius exceeds a critical value. In that regime, both the degree of asymmetry and 
the average expansion of the shock induced by the spiral mode increase monotonously 
with this ratio, exceeding the values obtained when a sloshing mode is artihcially im¬ 
posed. With a timescale of 2-3 SASI oscillations, the dynamics of SASI takes place fast 
enough to affect the spin of the neutron star before the explosion. The spin periods de¬ 
duced from the simulations are compared favorably to analytical estimates evaluated 
from the measured saturation amplitude of the SASI wave. Despite the simplicity of 
our setup, numerical simulations revealed unexpected stochastic variations, including 
a reversal of the direction of rotation of the shock. Our results show that the spin up 
of neutron stars by SASI spiral modes is a viable mechanism even though it is not 
systematic. 

Key words: hydrodynamics - instabilities - shock waves - stars: neutron - stars: 
rotation -- supernovae: general 


1 INTRODUCTION 

Neutron star spin at birth is a key parameter to associate 
pulsars and their progenitors. It carries information about 
the massive stars which explode in a core-collapse super¬ 
nova and give birth to a neutron star. Natal spins can be 
estimated via an extrapolation of the spin-down of observed 
pulsars and the determination of their current age. This 
indicates that a significant fraction of neutron stars are 
born with modest rotation periods in a broad range from 
a few tens to a few hundr eds milliseconds dPopov fc Turolla 
l2012l : lNoutsos et al.ll20lj ). Population synthesis studies sug¬ 
gest that the initial pulsar spin distribution may be mod- 
eled by a Gaussian distributio n centered around 300 ms 
dFaucher-Giguere fc Kas^l2006l) . However, early stellar evo¬ 
lution models found that much faster rotation close to 
breakup would result from the conservation of angular 
mom entum during core collapse (iHeger. Langer fc WoosleJ 
I2OOOI) . This large discrepancy in the distribution of pulsar 
spins at birth has to be resolved by better stellar evolution 


models and/or a better description of the role played by core 
collapse dynamics. 

Angular momentum transport processes throughout 
the stellar evolution play an important role by setting 
the rotation period of the core prior to colla pse. Recent 
aster o seismic observation s of red-giant stars (iBeck et al.l 
I2OI2I: iMosser et all l2012l) seem to require a more effi- 
cient angular momentum transport than us ually assumed 
(iDeheuvels et ah] l2014l : ICantiello et al.l l2014l ) . Considering 
the transport of angular momen tum by magnetic torques 
driven by a dynamo mechanism, iHeger. Wooslev fc SpruitI 
ll2005l) estimated pulsar spins of 10 — 15 ms - slower than 
previous estimates but still rather fast compared to observa¬ 
tional constraints. Transport of angular momentum by in ¬ 
ternal gravity waves ( see e. g. iTalon fc Charbonnell ll2003l) : 
iLee. Neiner &: Mathis! (l2014l) ) is also a p ossibility to explain 
the distribution of pulsar spins at birth. [Fuller et al.l (l2015l) 
estimated periods of 20 ms to 400 ms due to an influx of 
internal gravity waves during the latest burning stages. An- 
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gular momentum redistribution by hydrodynamic instabil¬ 
ities during the collapse has not been considered in these 
previous works. 

Detailed numerical simulations of core-collapse su¬ 
pernovae show that ma ssive stars do not ex plode 
in spherical symmetry ijLiebendorfer et al.l 1200 il l ex¬ 
cept for the low-mass end of sup ernova progeni¬ 
tors (IKitaura. Janka fc Hillebranda 200^. Neutrino - driven 
convection ( Herant et al.l 19941 : I Janka fc Muelleil Il996l l 
and the Standing Accretion Shock I nstability (SASI) 
l|Blondin. Mezzacaopa fc DeMarindl200.l| l are able to gener¬ 
ate large scale asymmetric motions and may play a crucial 
role in the success of the explosion (see e.g. iFoglizzo et al.l 
ll2015h for a recent review). Such asymmetric motions have 
been confirmed as the cons equence of a linear i nstability us¬ 
ing perturbatiye methods dFoglizzo et al.l[2007^ . The result¬ 
ing asymmetric explosions are supported by a series of ob- 
seryational eyid ences. The typical pulsar yelocities of a f ew 
hundreds km/s dArzoumanianrChernoff fc Cordeal2002l l is 
much higher than its progenitor ones. A pulsar kick imparted 
on the neutron star by a global deformation I = 1 has the 
potential to explain this yelocity distribution dScheck et al.l 
|2004 I2006I1 . Also, the distribution of '‘'‘Ti obseryed in 
Cassiopea A suggests a large scale asymmetric explosion 
dCrefenstette et al.ll2014h . 

Unstable modes of SASI can develop sloshing motions 
of the shock along a symmetry axis and also spiral motions 
when the axisymmetric constraint is released. SASI spiral 
modes can redistribute angular momentum during the col¬ 
lapse and significantly modify the neutron star spin from 
what could be inferred by ang ular momentum conservation 
dBlondin fc Mezzacappall2007fl . Using a simplified adiabatic 
model Blondin fc Mezzacappal ((2003) showed that a spiral 
mode can dominate the dynamics and is able to spin-up a 
neutron star to periods as short as 50 ms even if the progen¬ 
itor does not rotate. In this idealized scenario, some angular 
momentum is expelled in the explosion whereas the oppo¬ 
site angular momentum is accreted onto the surface of the 
neutron star. 

This mechanism relies on the dominant action of a 
spiral mode. In the linear regime of SASI a sloshing 
mode can be decomposed as two counter-rotating spi¬ 
ral modes with similar amplitudes and identical growth 
rates if the progenitor is non-rotating. A robust spiral 
mode may dominate the dynamics only if the symmetry 
between these counter-rotating spiral m odes breaks non- 
linear ly. Spira l modes we r e obt ained by iBlondin fc Shav^ 
(1200711 in 2D, iFernandezI 1I2OIOII in 3D using an approxi¬ 
mation of neutrino cooling. This numerical result has been 
confirmed by an expe rimental shallow-water analog of SASI 
(iFoglizzo et al.ll2012l l. Spiral modes dominate occasionally 
the dynamics in 3D numerical simulation of the collapse 
of a 27 Mq proge nitor using variou s approximations of 
neutr i no transp ort (|Hanke et a l1l2013l : ICouch fc O’Connoil 
I2OI4I : lAbdikamalov et al.l 12015117 Em ploying an analytical 
approach Guilet &: Fernanded (l2014h estimated the amount 
of angular momentum deposited when a single spiral mode 
dominates the dynamics in a non-rotating progenitor and 
showed that SASI has the potential to explain initial pulsar 
periods of a few tens to a few hundreds milliseconds. The 
slow end of this range of periods is compatible with the range 
100 ms -8s obtained bv lWongwathanarat. Janka fc Miillerl 


1 I 2 OI 3 II in their simulations of 15 Mq and 20 Mq progen¬ 
itors. An efficient spin-up mechanism driven by the SASI 
would requir e a long-lasting SASI activity up to the point 
of explosion. iMiiller. Janka fc Hegeri (I2OI2I I observed such a 
dynamics in their general relativistic neutrino hydrodynam¬ 
ics axisymmetric simulation of a 27 Mq progenitor. 

However, the non-linear dynamics of SASI may not al¬ 
ways be dominated by a spiral mode. Inde ed, no separa¬ 
tion o f angular momentum was obtained bv llwakami et aP 
(I2OO8II considering a model with both neutrino heating and 
cooling. Investigating the flow p a ttern below the shock wave, 
llwakami. Nagakura fc Yamadal (l2014l l showed that either 
sloshing or spiral modes dominate the dynamics depending 
on the mass accretion rate and neutrino luminosity consid¬ 
ered, illustrating the stochasticity of the angular momentum 
distribution in hydrodynamics simulations. Spin-up by the 
SASI can be effective only if the symmetry breaking leading 
to a single spiral mode has occurred before the explosion 
takes place. Even in this case, the amount of angular mo- 
mentu m accreted is sensitive to the position of the mass cut 
radius (iRantsiou et al.|[2011^ . 

We propose to study the timescale and rotational con¬ 
sequences of the symmetry breaking, which determines the 
respective roles of sloshing and spiral modes. We perform a 
set of 2D simulations of a simplified model of an accretion 
flow restricted to the equatorial plane, using cylindrical ge¬ 
ometry. The set of parameters considered shed some light 
on the non-systematic and non-deterministic features of the 
symmetry breaking. 

The paper is organized as follows. The physical and nu¬ 
merical models are described in section [2] Section |3] focuses 
on the properties of the symmetry breaking and the non¬ 
linear evolution of SASI to evaluate their consequences on 
the distribution of angular momentum. Our simulations are 
confronted to the dynamics of less idealized environments in 
section [T] in order to discuss the potential role of SASI on 
the initial neutron star spin. 


2 METHODS 
2.1 Physical model 

Our model consists of a standing accretion shock centered 
around a proto-neutron star in a stationary and non-rotating 
flow. For the sake of simplicity we focus our study on the 
equatorial plane of the collaps ing core, using cylindr i cal co - 
ordinates in a setup similar to I Yamasaki fc Foglizzo! 1I2OO8II . 
The main advantage of this model is to allow for non- 
axisymmetric modes of SASI in 2D. The accreting matter 
is modeled by a perfect gas with adiabatic index 7 = 4/3. 
Above the shock, the supersonic matter falls inwards radi¬ 
ally and reaches the shock radius rgh with an incident Mach 
number Mi = 5. Below the shock, the matter accretes sub- 
sonically onto the surface of the proto-neutron star which ra¬ 
dius is noted r*. A cooling function is included to mimic the 
neutrino emission due to electron capture with the approx¬ 
imation Co oc llBlondin fc Mezzacappal l2006fl . where 

p and P are respectively the density and the pressure. Neu¬ 
trino heating is neglected in order to suppress buoyancy in¬ 
duced convective motions and concentrate on SASI in its 
simplest form. The gravity is Newtonian and self-gravity is 
neglected. 
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The initial solution is computed by solving the time 
independent continuity, Euler and entropy equations be¬ 
low and above the shock. The two solutions are connected 
by the Rankine-Hugoniot jump conditions neglecting the 
dissociation of nuclei at the shock. The resulting dynam¬ 
ics only depends on the ratio of the initial shock to the 
proto-neutron star ra dii. Typical values of the radii ratio 
R = Tsh/r, are R 2 (ICouch fc O’Connorj|20l4l and i? « 4 
(iMarek fc Janka|l2009l ~l. depending on the progenitor struc- 
ture.Q When converting to physical units, we choose a proto¬ 
neutron star with radius r* = 50km and mass M* = 1.3Mq, 
and a constant mass accretion rate M = O.SMq s“^ as typ¬ 
ical values for the stalled shock phase of a core-collapse su¬ 
pernova during the first second after bounce. 


2.2 Numerical model 

To run our two-dimensional time-dependent hydro dynamic 
simu la tions we use a version of the RAMSES c ode (iTevssieil 
I2OO2I : iFromang. Hennebelle fc Tevssieil l2006fl adapted to 
cylindrical coordinates (r, <j>) for which the grid is uniformly 
spaced. RAMSES is a second-order finite volume code, which 
uses the MUSCL-Hancok scheme. The s imulations are per¬ 
form ed using the HELD Riemann solver (iMivoshi fc Kusand 
|2005^ and the monotonized central slope limiter. We set 
periodic boundary conditions in the azimuthal direction at 
(/) = 0 and 0 = 27r. The radial domain covers the interval 
r* ^ r ^ Tout with an outer boundary Vout/rsh ~ 3—6 so that 
the shock wave does not reach the edge of the domain. We 
impose reflexive inner bou ndary conditions and free oute r 
boundary conditions as in iBlondin fc Mezzacapp^ ll2006li : 
iFernandez fc ThompsonI ll2009bl ~l. Resolution effects are min¬ 
imized by fixing the number of radial cells below the initial 
shock to 150. The total number of radial cells is then in the 
range [540,1300] depending on the simulation. 1000 cells are 
used in the azimuthal direction, which is significantly larger 
than what can be afforded by current 3D simula tions (e.g. 
176 cells in iHanke et al.l[2013l : lMelson et al.ll2015fl . High res¬ 
olution is required to properly resolve the steep gradients of 
the flow dynamics in the vicinity of the proto-neutron star. 

An entropy cutoff is applied to the cooling function in 
order to avoid the divergence of the numerical solution in the 
vicinity of the proto-ne utron star ilFernandez fc ThompsonI 
l2009bl : lFernand^l2015l i. The cooling function is written 


C = Co exp 



( 1 ) 


where S = (7 — 1)~^ In {P/p'^) defines the entropy, Smin its 
value at r = r* and the value of k is chosen to introduce 
only minimal modifications to the stationary state flow. 

Once the numerical solution has relaxed on the grid for 
a few hundreds numerical timesteps, two density perturba¬ 
tions at pressure equilibrium are introduced ahead of the 
shock to trigger two counter-rotating spiral modes m = ±1 
ox m = ±2 as described in appendix]^ Perturbations are 


^ In these works, the ratio R is estimated by averaging the shock 
and neutrinosphere radii during the SASI domination phase of 
the dynamics. 


decomposed as Fourier modes in the azimuthal direction ac¬ 
cording to the general form 

M(r, t) = Re c^(r) j (2) 

where Cm{r,t) is the complex amplitude. We define the ini¬ 
tial e—asymmetry between spiral modes as: 

_ |cm| — |c-m| 

“|C^P + |C_^|2 

with |e| Si 1. Note that e = 0 corresponds to a mirror- 
symmetric sloshing mode and e = ±1 to a single spiral mode. 

Our aim is to estimate the timescale for a symmetry 
breaking, after the phase of linear growth which lasts less 
than ~ 3 SASI oscillations (< 100 ms). Two different meth¬ 
ods have been developed for that purpose. The first one is 
based on the time evolution of the angular momentum flux 
through the inner boundary. This flux is very close to zero 
for a sloshing mode and starts to deviate from zero once one 
of the spiral modes dominates. The second method is based 
on the angular tracking of the minimum shock radius. This 
point corresponds to one of the triple points that form in 
the shock wave. Its rotation rate evolves rather erratically 
for a sloshing mode but becomes fairly constant for a spiral 
mode. The two methods are consistent within a SASI period 
which is sufficient for our study. 

Our code has been carefully tested to check that it does 
not introduce any artificial source of asymmetry. If the ini¬ 
tial density perturbation is mirror-symmetric, i.e. e = 0 , 
the mirror-symmetry is conserved and the sloshing mode 
oscillates along a fixed axis. Moreover, two simulations with 
opposite initial perturbations, e and —e, show two dynami¬ 
cal evolutions that remain mirror-symmetric within machine 
precision. Besides, the robustness of the code has been tested 
by comparing the growth rates and oscillatory frequencies 
measured in our s imulations to those obtained with a pertur¬ 
bative analysis bv lYamasaki fc Foglizzd (l2008l j. The discrep¬ 
ancies are less than 8 % for the growth rates and less than 
2% for the oscillatory fr equencies. This is similar t o the g ood 
agreement obtained by [Fernandez fc ThompsonI (l2009bf) . 


3 RESULTS 

3.1 A critical ratio for the symmetry breaking 

We performed a total of 80 simulations vary¬ 
ing the two parameters R and e such that R = 
{1.67, 2, 2.22, 2.5, 3, 3.5, 4} and 10"® ^ je] ^ 1. Our 
simulations sho w that contrary to wha t was o btained in 
some studies llBlondin fc Mezzacappai l2007l : iFernand^ 
l 2010 ll . a spiral mode does not always dominate the late 
evolution. The ratio R was found to determine whether the 
symmetry breaking occurs or not. When R ^ 2, the late 
evolution is dominated by a robust sloshing mode, even if 
a single spiral mode (i.e. |e| = l) was used to perturb the 
stationary flow (Figure [T] left). The azimuthal index of the 
sloshing mode can either be m = 1 or m = 2 , depending 
on the value R. In this regime, angular momentum is not 
significantly redistributed. 

A totally different behavior is observed when R > 2. 
A spiral mode dominates the late evolution (Figure[T]right) 
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Entropy at t = 960 ms Entropy at t = 187 ms 



-200 -100 0 100 200 -200 -100 0 100 200 
X (km) X (km) 


Figure 1. Left: Entropy snapshot at t = 960 ms for R = 2 and e = f. A sloshing motion dominates the non-linear regime despite a spiral 
perturbation. Right: Entropy snapshot at t = 187 ms for R = 3 and e = 0.1. The symmetry breaking has already occurred. (Animated 
versions of these figures are available in the online journal.) 


even for weak e—asymmetry, enabling a redistribution of 
angular momentum. These results raise the question of the 
mechanism responsible for this symmetry breaking. The dy¬ 
namics of SASI observed in our simulations may help to 
characterize this mechanism as discussed in section im 


3.2 Timescale for the symmetry breaking 

We apply the methods described in section 12.21 to compute 
the timescale to reach a symmetry breaking as a function 
of i? > 2 and e (Figure (H). For R = {2.5, 3, 4} the symme¬ 
try breaking occurs within 2 to 10 SASI oscillations in the 
non-linear phase, which is fast enough to potentially redis¬ 
tribute angular momentum before the explosion. In the case 
R = 2.22, a spiral mode dominates only after 20 to 30 SASI 
oscillations. This timescale may be too slow to impact the 
neutron star spin. The ratio R = 2.22 illustrates the conti¬ 
nuity between a rapid symmetry breaking and an absence 
of symmetry breaking. 

The influence of the initial asymmetry on the timescale 
is not straightforward. If the asymmetry is large enough (i.e. 
|e| ^ 0.2), the timescale decreases with |e| as would be in¬ 
tuitively expected. The trend seems rather chaotic and the 
uncertainties on the timescale are as large as the variability 
of the results when |e| ^ 0.2. Furthermore, the direction of 
rotation of the spiral mode is not always the one determined 
by the initial asymmetry. Indeed, approximately half of our 
simulations with |e| ^ 0.2 show a symmetry breaking in 
the other direction (square symbols in Figure (2)1. The code 
has been extensively tested to prevent numerical artifacts 
from inducing asymmetries. This non-deterministic feature 
could instead be generated by several non-linear processes 
which we mention as possible paths towards an explana¬ 
tion. The first one is based on the parasitic instabilities, 
such as Kelvin-Helmholtz and Rayleigh-Taylor that have 
been proposed to explain the satur ation amplitude of the 
SASI (iGuilet, Sato fc Foglizzo! l2010l ). The parasites which 
develop on SASI spiral modes might modify the asymmetry 



Figure 2. Number of SASI oscillations before reaching a symme¬ 
try breaking with respect to the initial asymmetry e. From top 
to bottom are shown ratios R = 2.22, 2.5, 3, 4. Square symbols 
show cases for which the direction of rotation is opposite to the 
one of the initial asymmetry. 


level in a stochastic way before the symmetry breaking (Fig¬ 
ure EJ. The second non-linear process relies on secondary 
shocks that arise before the symmetry breaking. Multiple 
secondary shocks, shown in F igure jd] were witnessed by 
iFernandez fc ThompsonI ll2009ah . These shocks may interact 
with the global advective-acoustic cycle. The entropy and 
the vorticity produced by secondary shocks may produce an 
acoustic feedback in the azimuthal direction opposite to the 
initial acoustic wave. This phenomenon might be able to al¬ 
ter the competition between counter-rotating spiral modes 
and add some stochasticity before a symmetry breaking oc¬ 
curs. An example is shown in Figure [S] where a secondary 
shock is able to generate opposite angular momentum well 
inside the outer shock wave. 
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Figure 5. Snapshots of the entropy (left) and the angular momentum (right) at t = 942 ms for R = 4 and e = 1 limited to the inner region 
of the domain (see Figure [4| for a broader view). The red lines display the location of the secondary shocks. A secondary shock generates 
higher entropy material. This region corresponds to positive angular momentum material (red regions below the secondary shock) whereas 
the dynamics is dominated by a spiral mode containing negative angular momentum (blue regions in the angular momentum snapshot). 



Figure 3. Entropy snapshot at £ = 112 ms for R = 3 and 
e = 0.01. Both Kelvin-Helmholtz (“KH”) and Rayleigh-Taylor 
(“RT”) structures are visible after 3 SASI oscillations. These in¬ 
stabilities may add stochasticity before the symmetry breaking 
between SASI spiral modes. 

3.3 Reversal of the direction of rotation 

This section and the following one are dedicated to measur¬ 
ing the rotation induced by the spiral mode. We focus on a 
set of 7 simulations where the radii ratio R is varied and the 
initial asymmetry is set to e = 1. The angular momentum 
density profile 

lz{r,t) = j ( 4 ) 

is averaged in time over ten SASI periods in the non-linear 
phase as shown in Figure [6] Two counter-rotating regions 
are observed and the radius separating them is labeled ro. 


Figure 4. Snapshot of |V (Pr^)| /Pr® at i = 942 ms for P = 4 
and e = 1. Multiple secondary shocks are present. The solid white 
lines delimit the domain zoomed in Figure [5] 

If there is no symmetry breaking (P ^ 2), the average pro¬ 
file is best defined using the linear phase only. The angular 
momentum Lz{t) contained between the radius ro and the 
shock wave is computed by: 

/■rsh 

Lz{t) = / lz{r,t)dr. (5) 

Jrg 

Figure [3 shows the time evolution of the enclosed angu¬ 
lar momentum for our set of simulations. The cases R = 4 
(solid red line) and R = 2.22 (dashed black line) exhibit a 
surprising inversion of the direction of rotation of the spiral 
wave. Both events take place in the fully non-linear regime. 

For R = 2.22, the change of direction lasts approxi¬ 
mately 8 SASI periods during which a sloshing mode dom¬ 
inates. The angular momentum produced by SASI is very 
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Figure 6. Time averaged angular momentum density profile in 
the non-linear regime for R = S normalized by 



t(s) 


Figure 7. Time evolution of the enclosed angular momentum 
normalized by for 7 different values of R. 


Figure 8. The analytical estimate of the initial neutron star spin 
period (dot-dashed blue line) is compared to the non-linear regime 
of our simulations (red bars). The bars refer to the time varia¬ 
tion of the amount of angular momentum accreted. For R < 2.5, 
various non-linear effects can cancel the angular momentum re¬ 
distribution and lead to very slowly rotating neutron stars. 

Figure [S] shows a comparison between the analytical esti¬ 
mate (Ea. lB19l using for the spiral mode amplitude the value 
measured in the simulation as described in next subsection) 
and the time averaged value in the non-linear regime of our 
simulations. Our results are consistent with the analytical 
estimates within a factor 2 for R > 2.5 and confirm that 
spiral modes of SASI are able to spin-up a neutron star to 
periods of tens to hundreds milliseconds. The larger discrep¬ 
ancies for i? < 2.5 are no surprise because the spiral modes 
do not exist in the non-linear regime [R = {1.67, 2}) or 
a reversal of the direction of rotation takes place during a 
significant fraction of the non-linear regime {R = 2.22). For 
that range of parameters, angular momentum redistribution 
by SASI is inefficient to spin-up the neutron star. 


low during that period. Even though a symmetry breaking 
has occurred, we observe that the non-linear dynamics of the 
SASI is able to cancel the angular momentum redistribution 
for a significant time. In the case i? = 4, the change of the 
direction is achieved on a much shorter timescale, less than 
a SASI period. 

The robustness of this behavior was confirmed by re¬ 
peating these two puzzling simulations varying slightly one 
parameter such as the numerical resolution or the pertur¬ 
bation amplitude. This intriguing phenomenon calls for a 
physical interpretation. A possibility might be that the sec¬ 
ondary shocks discussed in Sect. 13.21 break the advective- 
acoustic cycle and establish temporarily a new one between 
the second shock and the feedback region. This process is il¬ 
lustrated by Figure [5] The conditions for this adverse contri¬ 
bution to be able to reverse the direction of rotation remain 
to be determined. 

3.4 Estimate of the pulsar spin 

iGuilet fc Fernandej (|2014h derived an analytical estimate 
of the angular momentum redistribution driven by a single 
spiral mode in spherical geometry. This approach has been 
adapted to the cylindrical geometry in appendix [B] 

Birth periods of neutron stars are inferred from our sim¬ 
ulations using a moment of inertia of 7 = Ins x 10'*®gcm^. 


3.5 Saturation amplitude of SASI 

The saturation amplitude of SASI Ar is a key element of 
the spin-up by spiral modes because the amount of angu¬ 
lar momentum redistributed scales as Ar^ (Eg. IBlSl) . The 
increase of the saturation amplitude with the ratio R (Fig¬ 
ure [§] left) is consistent with the highest spin obtained for 
highest values of R (Figure [H|. 

However, the satu ration amplitude obtained by ap¬ 
plying the formalism of iGuilet. Sato fc Foglizz 3 ll 2010 l) de¬ 
creases with increasing R (Figure [5] left). The higher satu¬ 
ration amplitudes observed at large values of R in the sim¬ 
ulations indicate that in this regime the parasitic instabil¬ 
ities are not as efficient at stopping the gro wth of SASI as 
predicted bv iGuilet. Sato fc Foglizzg (l2010l) . This suggests 
either that a more elaborate description of the parasitic in¬ 
stabilities is necessary or that another process is responsible 
for the saturation of SASI. 

The shock expansion due to SASI is increasing with 
R more steeply than the saturation amplitude: between 
R = 2 and R = 4, it increases by a factor ~ 4 
while Ar/rgh increases by a factor ~ 2 (Figure O. This 
is consistent with the shock expansion varying quadrati- 
cally with the saturation amplitude as might be expected 
for a non-linear effect. A simil ar trend was observed by 
iFernandez fc ThompsonI ll2009bl ) with a slighter increase 
which may be attributed to the geometry difference. A di- 
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Saturation amplitude of SASI 


Shook expansion 




R 


R 


Figure 9. Left: The satur ation amplitude of SASI c ompu ted by 
applying the formalism of iGuilet, Sato &: Foglizza 1I2OICIII (dot- 
dashed blue line) is compared to the one in our simulation for 
e = 1 and 7 different values of R. In the simulations the satu¬ 
ration amplitude is estimated by averaging the amplitude of the 
dominant mode over the non-linear regime and normalized by 
the initial shock radius (solid black line) or by the average shock 
radius in the non-linear regime (dashed black line). Right : Vari¬ 
ation of the shock radius compared to its initial value as a 
function of R for e = 1. 


rect comparison with si mulations of steady-state flows in¬ 
cluding neutrino hea ting JOhnishi, Kotake fc Yamada||200fil: 
llwakami et al.l2008l : llwakami. Nagakura fc Yamadal 20141 ) is 

less straightforward because R is also affected by the neu¬ 
trino luminosity. Larger ratios may correspond to dynamical 
evolutions dominated by neutrino driven convection. 


3.6 Differences between spiral and sloshing modes 


For R > 2 the saturation properties of the spiral mode 
(e = 1) are compared with the ones of a sloshing mode ob¬ 
tained by imposing mirror-symmetric initial perturbations 
(e = 0) (Figure fTOl left). If the ratio R is close to the thresh¬ 
old for symmetry breaking, the average shock radius and the 
saturation amplitude are almost equal between a sloshing 
mode and a spiral mode. When the ratio R is large enough 
for the domination of a spiral mode, the shock radius and 
the saturation amplitude are increased by up to about 40% 

compared to the mirror symmetric evolu tion. _ _ 

These results confirm the work of iFernand^ (l2015l ) 
which showed that a spiral mode in 3D may lower the criti¬ 
cal neutrino luminosity compared to a sloshing mode in an 
axisymmetric case that generates less non-radial kinetic en¬ 
ergy. In our simulations with the highest ratios R, spiral 
modes are indeed able to double the total non-radial ki¬ 
netic energy compared to sloshing modes (Figure [Till right). 
as observe d in th e simulations without neutrino heating of 
[Fernand^ ll2015h . For smaller ratios R on the other hand, 
the difference of kinetic energy between spiral and sloshing 
modes is more modest, suggesting the existence of differ¬ 
ent regime s. This might be t he reason why other groups, in 
contrast to [Fernand^ ll2015l) , had not found a lowered crit- 
ical luminosity in 3D simulations exhibiting a spiral mode 
dHanke et al.l 20131 ). 


Saturation amplitude, shock expansion Non-radial kinetic energy 




Figure 10. Left: Ratio of the saturation amplitude (dashed black 
line) and the mean shock radius (solid blue line) between a spiral 
mode and a sloshing mode for 5 different values of R. Right : Ratio 
of the average non-radial kinetic energy in the non-linear regime 
between a spiral mode and a sloshing mode. Note that spiral 
modes are systematic outcomes for the range of R considered 
in this figure, unless a mirror-symmetric sloshing is enforced by 
choosing e = 0. 


3.7 A possible path towards the symmetry 
breaking mechanism 

A description of the physical processes responsible for the 
symmetry breaking would be helpful to anticipate the effi¬ 
ciency of SASI at spinning-up neutron stars in more realistic 
models. A first constraint is that no spiral mode dominates 
the non-linear dynamics if i? ^ 2. Additional clues may be 
inferred from the following properties: 

• Unlike for m = 1 modes, the m = 2 sloshing mode 
in our setup is never transformed non-linearly into a spiral 
mode. 

• In the case R = 2, despite a linear domination of the 
mode m = 2, the mode m = 1 eventually prevails due to a 
non-linear coupling between these modes (Figure [TT}. How¬ 
ever, the transition between these modes does not lead to a 
spiral mode (Figure [T] left). 

• Linearly, the mode m = 2 dominates the mode m = 1 
for R ^ 2.2. Interestingly, the critical ratio for the symmetry 
breaking is close to this linear transition. 

• The efficiency of the symmetry breaking seems to be 
linked to the difference of saturation amplitudes between 
the spiral and the sloshing modes. A fast symmetry breaking 
corresponds to a significantly larger amplitude of the spiral 
mode, while the amplitudes are approximately equal when 
the symmetry breaking is slow {R — 2.22, see Figures [2] and 
[To] left). 


4 DISCUSSION 

Several simplifications have been made in our model to study 
the physics of SASI in its simplest form and less idealized 
models might modify some aspects of our results. The di¬ 
mensionality and the geometry are important points to raise. 
The density and velocity profiles in cylindrical geometry are 
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t (s) 


Figure 11. Time evolution of the amplitudes of the Fourier 
modes m = 1 (thin blue line) and m = 2 (thick green line) 



r (km) r (km) 


Figure 12. Density profile (left panel) and radial velocity profile 
(right panel) computed in cylindrical coordinates (solid red line) 
and in spherical coordinates (dashed blue line). 


compared in Figure[T2]to those obtained in spherical geom¬ 
etry for the same parameters used in Sect. [2l In the sub¬ 
sonic region of the flow, the density profile is independent of 
the geometry but the advection time is shorter in spherical 
geometry. Remembering that SASI frequencies and growth 
rates scale like the advection rate, these quantities are higher 
in spherical geometry. A symmetry breaking between SASI 
spiral modes may therefore occur earlier in a 3D spherical 
model than in 2D cylindrical geometry. The dimensional¬ 
ity of the model impacts the amount of angular momentum 
via its dependence on the saturation amplitude. The lat- 
ter was found to be weakly sensitive to the dimensi onality 
llFernandejl2O10| : lHanke et al.|[201^ : lFernandejl2015l ). How¬ 
ever, drawing conclusions on this issue may require to clarify 
the divergence between the predicted and the measured sat¬ 
uration amplitudes observed in our study (Figure [H]). 

The initial rotation of the progenitor has been ne¬ 
glected for the sake of simplicity, but could dominate the 
angular momentum budget if it is fast enough. Considering 


the development of spi r al mo des in a rotating progenitor, 
iBlondin Sz Mezzacappal (l2007l ) showed that SASI could sur¬ 
prisingly decelerate a neutron star which accretes SASI in¬ 
duced angular momentum opposed to the initial rotation of 
the stellar core. They also showed that this mechanism could 
even lead to the forma t ion of a counter rotating neutron star. 
lYamasaki fc Foglizzol (l2008l ) confirmed that rotation favors 
prograde spiral modes and showed that SASI growth rates 
depend linearly on the angular momentum of the progeni¬ 
tor. These results raise the issue of the critical rotation rate 
of the progenitor above which the neutron star spin at birth 
is mostly determined by the conservation of initial angular 
momentum. A crude estimate can be made by evaluating 
the angular momentum accreted by a neutron star during 
the collapse of a non rotating core. However, the effect of 
rotation on the saturation amplitude of SASI is still poorly 
known. The amount of angular momentum redistributed by 
a spiral mode may increase even if the centrifugal force is 
negligible. The mutual influence of the initial rotation and 
the SASI induced dynamics on the birth period of neutron 
stars will be addressed in a forthcoming paper. 

Taking into account neutrino heating would add a 
source of stochasticity through the development of neutrino- 
driven convect ion, and help address the diversity of ex¬ 
plosio n paths (iFernandez et al]l2014l : ICardall fc Budiardial 
l2015h . Pre-collapse convective asymmetries may also add 
stochasticity to post-bounce dynamics and affect the fate o f 
the massive star (ICouch fc C)tll20ld : Miiher fc JankalboiSl ) . 

Ilwakami. Nagakura fc Yamadal (2014) explored the di¬ 
versity of flow patterns behind the stalled shock and ob¬ 
served that the symmetry breaking is not systematic in 
their SASI-dominated mode l A with M — 0.2 Mf7|S~^, Li, = 
2x lO^ ^erg s~^ (see table 1 in flwakami. Nagakura fc Yamadal 
(|2014|) ). Changing slightly the numerical resolution or the 
noise in the initial conditions either lead to a quasi¬ 
stationary sloshing or spiral mode. However, their other 
SASI dominated models exhibit only spiral modes. This 
opens the question of the existence and the value of a critical 
ratio R for the symmetry breaking in more complex models 
which may be more subject to stochasticity. Wit hout neu¬ 
trino heating in 3D simulations iFernand^ (l20ld) observed 
that the amount of angular momentum redistributed by the 
dynamics of SASI is greatly reduced for R = 1.67 compared 
to R = 2. These results are consistent with the fact that 
J? = 1.8 is the transition betw een I = 1 and I = 2 lin¬ 
ear modes llFoglizzo et al.|[2007^ . A similar transition takes 
place in cylindrical geometry for Rk. 2.2. 


5 CONCLUSION 

A simplified setup in cylindrical geometry has been used 
to investigate the flow pattern in the non-linear regime of 
SASI for a non-rotating progenitor. A symmetry breaking 
between counter rotating spiral modes occurs only if the ra¬ 
tio of the initial shock to neutron star radii R > 2. If this 
condition is satisfied, the dynamics is dominated by a spiral 
mode, independently of the initial conditions and SASI has 
the potential to spin-up a neutron star to initial periods of a 
few tens to a few hundreds milliseconds. However, if R 2, 
there is no sign of symmetry breaking and a sloshing mode 
dominates the dynamics. This case leads to very slowly ro- 
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Spin-up by SASI spiral modes 9 


tating neutron stars (Figure [8]). These properties set strong 
constraints on the still unknown mechanism responsible for 
the non linear symmetry breaking. 

The timescale for symmetry breaking of the order of 
2-3 SASI oscillations is short enough to affect the angular 
momentum of the neutron star before the explosion. This 
timescale shows stochastic variations when the initial asym¬ 
metry is weak (Figure [2j. Memory of the initial perturba¬ 
tions is lost before a symmetry breaking can occur. Moreover 
the non-linear dynamics of the SASI can lead to a change of 
the direction of rotation (Figure [7|. This unexpected result 
reveals how complex the dynamics can be even in a simpli¬ 
fied setup. Additional sources of stochasticity are expected 
from the development of neutrino-driven instabilities as well 
as pre-collapse convective asymmetries. 

Spin-up by SASI may lead to birth periods of 
neutron s tars compatible with obse r vations as pro¬ 
posed by llBlond in fc M ezzacaop^ l2007l : iFernand^ I 2 OI 0 I : 
lOuilet fc Fernandej l2014l 'l. The neutron star periods ob¬ 
tained in our simulations are cons istent with analytical 
estimates llGuilet fc Fernandej l20f3 ~) regardless of the di¬ 
mensionality of the setup considered. Diverging conclusions 
regarding the efficie n cy of the spin-up me chanism (e.g. 
Ilwakami et al.l (l2008l i: iRantsiou et al.l (l201ll li may be ex¬ 
plained by the choice of progenitors or parameters favoring a 
dynamics dominated by neutrino-driven convection instead 
of SASI. 

If the shock radius is large compared to the neutron 
star radius, the saturation amplitude of the spiral mode 
is larger than the sloshing mode resulting from a mirror- 
symmetric evolution. On the one hand, the larger kinetic 
energy and shock expansion induced by spiral modes would 
presumably support a lo wered critica l neut rino luminosity 
in 3D, in agreement with iFernande j (l2015l ). On the other 
hand, when the shock is closer to the neutron star, the sat¬ 
uration amplitude of the spiral mode becomes closer to that 
of the sloshing mode IFigure fTCTl) . In that regime, the ax- 
isymmetric and 3D dynamics may be expected to be more 
similar. 

The highest saturation amplitudes observed for a large 
shock radius seem har d ly ex plained by the formalism of 
lOuilet, Sato fc Foglizzol (l2010l ) and require further investi¬ 
gations. 

Initial rotation in the stellar core has been neglected in 
this study in order to focus on the rotation induced by the 
spiral mode of SASI. Rotation is more likely to trigger spiral 
mode s llBlondin fc Mezzacaopal l2007l : lYamasaki fc Foglizzol 
l2008l 'l which may spin-down the neutron star or even give 
birth to a counter rotating one. Nevertheless, the impact 
of the rotation on the saturation amplitude of SASI and 
its non-linear dynamics is still poorly known. Such a study 
would help characterize how SASI can affect the mapping 
between the angular momentum profile of massive stars and 
the distribution of the initial pulsar spins. The influence of 
initial rotation in the core will be addressed in a forthcoming 
paper in order to disentangle the respective contributions of 
the initial angular momentum and the dynamical effects of 
SASI on the pulsar spin at birth. 
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APPENDIX A: DENSITY PERTURBATIONS 

In the linear regime, the two spiral modes of index ±m are 
triggered by over-densities injected at the outer boundary of 
the domain: 

5p± {6, t) = A{1 ± e) cos {tUr t^m6) (Al) 

where A and cur are the amplitude of the perturbations and 
the oscillatory frequency of the SASI mode, respectively. In 
this formulation — 1 ^ e ^ 1 and the sign of e selects the 
dominant spiral mode in the linear regime. The overall per¬ 
turbation are written as: 

6p {6, t) = H{t) {5p+ (9, t) + Sp. {9, t)) (A2) 

where H(t) is a function used to smoothen the perturbation 
such that: 

exp |- j^ ^-fa+jadv/4) j^| if ^ i ^ 

4 ^ ^ 5 7”adv/4 

t-(to-t5r,dv/4) j 'I if 5r,dv/4 < t ^ 6r,dv/4 
,0 otherwise 

(A3) 

where to is the time when the perturbations start to be ad- 
vected through the outer boundary, Tadv = 2 7r/o;r is the 
advection time and a is a coefficient used to vary the ampli¬ 
tude of H{t) from 10“^® to 1 over a timescale 


H{t) = 


1 if "^adv/ 

expj- 


APPENDIX B: ANGULAR MOMENTUM 
REDISTRIBUTION BY SASI SPIRAL MODES 
IN GYLINDRICAL GEOMETRY 


In this appendix, we ada pt the formalism developed by 
iGuilet fc Fernandej ll2014ll for the angular momentum re¬ 
distribution by a SASI spiral mode in spherical geometry, to 
the cylindrical setup considered in this paper. Because the 
derivation of the equations follows very similar steps, we do 
not reproduce all of them here but highlight the differences 
linked to the change of geometry. The end result takes a very 
similar form to the spherical geometry, differing only in the 
numerical factor. 

As in the rest of the paper, we consider a 2D accre¬ 
tion flow in cylindrical geometry, assumed to be invariant 
in the vertical direction and described using cylindrical co¬ 
ordinates {r, (f)}. The surface integrated angular momentum 
density is defined in equation Q (where p is to be under- 
stood as a vertically integ rated surface density). While in 
iGuilet fc Fernandez (12014 1 the surface integration was done 
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on spherical shells, it is here performed on cylinders. An- 
Rular momentum conserva tion can then be written as in 
lOuilet fc Fernandej (|2ni4l ~l 


Defining TiieyO,m as the Reynolds stress amplitude of a given 
mode with Fourier index m and with the time dependence 
scaled out, 


dth + drT = 0 , 


(Bl) 


where J- is the angular momentum flux integrated over a 
cylindrical surface 

T{r, t) =F f pVrV^ dcj). (B2) 


As in iGuilet fc Fernandej ll2014l ~). the flow is described 
as a stationary background with superimposed small ampli¬ 
tude perturbations 


p{r,(j),t) 

= po{r) + &p{r,4>,t) + &'^p{r,4>,t)A... 

(B3) 

Vr 

= Wo + <5fr + + ... 

(B4) 

Vcjy 

= Sv^ + S'^V^ + ■■■ 

(B5) 


where <5 and denote first- and second order Eulerian per¬ 
turbations, respectively, with 5^5^. With this decomposi¬ 
tion, the surface integrated angular momentum density and 
flux read 


4 


Mr r Sp Sv^ ^ 5‘^v^ 
‘2-n- J [po Vo Vo 
IzVo + Tfley 


d(j), 


(B6) 

(B7) 


where M = —2nrpoVo is the stationary mass flux, and Tji^y 
is the surface-integrated Reynolds stress 


rr ! f &Vr&V^ , , 

TReyir,t) = --- / - 2 ^d(l>. (B8) 

First order perturbations are then decomposed into a super¬ 
position of spiral modes with sinusoidal angular dependence 
with Fourier index m fthis replaces the spherical harmonics 
decomposition used in I Guilet fc Fernandej (l2014h 'l. and the 
time-dependence of a plane wave with complex frequency 
u! = oJr + ivJi, with Ldr and iVi the real and imaginary parts, 
respectively. The space and time dependence of an arbitrary 
first-order perturbation 5A is therefore 

5A{r, (j), t) = ^ Re j (B9) 


where 5Am{r) is the complex amplitud(j3- The radial struc¬ 
ture and eigenfrequencies of these modes can be comp uted 
with a linear analysis as in I Yamasaki fc Foglizzol (l2008ll . 

Using this linear eigenmodes decomposition, the 
Reynolds stress can be written 


Tuey 


Mr Vo 
2 




Sv4, ,7n ^^r,Tn 

Vo Vo 



(BIO) 


Mrvo 

TneyQ,7n\V) — - — - Rc 


^V(j)^ni, ^^r,77 

Vo 


we can write equation (|B10|) as 

TRey — 


(Bll) 


(B12) 


Gombining equations JbU, (Ezl), and (IB12|l . angular 
momentum conservation is written as a partial differential 
equation for the evolution of Iz 


dtlz + drihvo) = -Y, drTR,yO,me^‘"'’^\ (B13) 

m 


The angular momentum density R can therefore be written 
as the sum of contributions from different Fourier mo des, 
with each term given by fsee lGuilet fc Fernandejl2014l . for 
more details on the derivation) 


^zO.m — 


Vo 






Vo 


rr 


rj-, I 

— -i-ReyO.m ur, 


^^0 


(B14) 


where radv(r') = ^ dr/vo is the advection time from the 

shock radius r^h to a radius r < rsh. 


Bl Angular momentum density below the shock 


The angular momentum density below the shock due to a 
spiral SASI mode follows from equations (IBllll and (IB14II . 


^zsh — 


TReyoGsh) _ Mrsh 


Re 


I 5V4,6Vr 




(B15) 


This expression can be evaluated using the boundary condi¬ 
tions of linear eigenmodes at a shock with a constant disso¬ 
ciation energy 


Mvsh 


^rVsh. 

27rUsh 


/(«,A4i) 



2 


(B16) 


where Ugh is the radial velocity below the shock, Ar is the 
amplitude of the shock deformation induced by the SASI 
spiral mode, ac is the compression ratio of the shock, A4i is 
the upstream Mach number, and /(«:, Af i) is a dimensionless 
factor 


= TT (k - 1) (1 - 1/k) 


1 + l/Ml 


7 — (7 -I- 1)/k + l/M\ ' 


(B17) 


^ Contrary to the spherical case, where the azimuthal velocity 
perturbation 6V(f) has a different angular dependence than the 
rest of the variables, the cylindrical geometry allows to describe 
all the va riables with the same Fourie r decomposition. When com¬ 
paring to lGuilet &: Fernandej ll2014lb this difference of decompo¬ 
sition leads to a numerical factor im being included into the com¬ 
plex amplitude of the azimuthal velocity. This —together with the 
different normalization of Fourier modes compared to spherical 
harmonics— is the reason why the Reynolds stress expressions in 
equations llBTiil and pSllJ differ by a fac t or im /A-n from equa¬ 
tions (17) and (18) of lGuIlet &: Fernandej ll2014h . 


B2 Approximate expression for the angular 
momentum contained in a spiral wave 

Using the same method as iGuilet fc Fernandej (l2014l ) we 
obtain an approximate expression for the total angular mo¬ 
mentum contained in the spiral wave (i.e. between the shock 
and the radius where the angular momentum changes sign) 

Lz ~ {rsh-rt.)lzsh ^ m/(K, — p^Mr^h^ (“) • 

27r|'t^sh| \t’sh/ 

(BI 8 ) 
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Spin-up by SASI spiral modes 11 


This is the same equation as iGuilet fc Fernanded ll2014l '). but 
note that the numerical factor differs by a factor 

47 r. This is mostly due to the different normalization of the 
Fourier modes co nsidered here compared to the spherical 
harmonics used in iGuilet fc Fernanded (l2014h . Considering 
the moment of inertia I of the neutron star, this can be 
translated into a minimum period of uniform rotation 

27r|ush| 1 (B19) 

mf{K,Mi) ujr{rsh-rt) Mrsh'^ 
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